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A highly extrapolatable semiempirical shell model mass equation applica- 
ble to translead elements up to Z = 126 is presented. The equation is applied 
to the recently discovered superheavy nuclei 293 118 and 289 114 and their decay 
products. 
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^ ; I. INTRODUCTION 

On- 

A recent experiment [|TJ is consistent with the formation of the nucleus 293 118 and its 

sequential decay down to 269 Sg (Z = 106). The a-decay energies vary rather smoothly along 

the chain, precluding the traditional macroscopic-microscopic [0,|3] Z = 114 as a major magic 

proton number in these nuclei. Recent phenomenological studies of B(E2) M and Wigner 

term || systematics indicate Z = 126 as a plausible next spherical proton magic number after 

lead. Recent self-consistent and relativistic mean field calculations [§pf| variously predict 

proton magicities for Z = 114, 120, 124 and 126, depending on the interaction used. 

Contrary to most of the above findings, the semiempirical shell model mass equation 
(SSME) U is based on the assumption that Z = 114 is the next proton magic number 
after lead, and it stops there. Moreover, the quality of its agreement with the data starts 
deteriorating already beyond Hs (Z = 108). (See fig. 1, and fig. 4 of ref. ||.) One has to 
find a substitute for the equation in the neighbourhood of Z = 114 and beyond. 

In the early stages of developing the SSME [TIJ both Z = 114 and Z = 126, then consid- 
ered possible alternative candidates for the postlead proton magic number, were tried as an 
upper shell-boundary for translead elements. The agreement with the data was about the 
same for both choices, and considering the prevailing view in the mid nineteen-seventies Z = 
/\ 114 was chosen for the SSME mass table. We have recently [|ll[j established a high predictive 



O 

o 



power of the early Z == 126 results in the interior of the shell region with Z > 82 and 126 
< N < 184 (called here region B) by comparing them to the newer data measured since 
then, and proposed using them as a substitute for the SSME M in superheavy elements 



(SHE) research. We have also [12] established a high predictive power for Q a values of the 
early Z = 126 results in the interior of the shell region with 82 < Z, N < 126 (called region 
A). However, the quality of the predicted masses and other mass differences worsened much 
compared to that of the original adjustment. Therefore we readjusted the coefficients which 
largely cancel in Q a , and proposed using the resulting equation as a substitute for the SSME 
II in the interior of region A. 



* Contribution to the 2nd Euroconference on Atomic Physics at Accelerators: Mass Spectrometry, 
Cargese, 19-23 September 2000. 

^Present address: Kashtan 3/3, Haifa 34984, Israel 



In sect. 2 we give the Z = 126 equation in the two regions and briefly discuss its predictive 
power. In sect. 3 we apply and comment on it in relation to SHE research in region B. 
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FIG. 1. Deviations of the mass predictions [9] from the data for the 56 new masses in region B 
measured after the original adjustments were made. 

II. THE MASS EQUATION AND ITS EXTRAPOLATABILITY 



In the SSME the total nuclear energy E is written @,0 as a sum of pairing, deformation 
and Coulomb energies: 



E (N, Z) = E pair {N, Z) + E def {N, Z) + E Coul (N, Z) 



(1) 



/2Zn\ 1/3 



Ecoui (N, Z) = \-^j [a c + f3 c (Z- Z ) + 1° (Z - Z f 



(2) 



E pair (N,Z) = ( r j)[« + /3{A- A Q ) + 1 (A- A f + eT (T + 1) 
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for region A, and 

E pair (N, Z) = (^j [a + (3 X (N- N ) + (3 2 {Z - Z ) 

+ 7l (7V - N ) 2 + 72 (Z - Z ) 2 + 73 (iV - N )(Z - Z ) 
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for region B, 
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E def (N,Z) = (-£) [ipn$u (N, Z) + ^ 20 [^ 20 (N, Z) + ^ 20 (Z, N)\\ 



(3) 



(4) 



(5) 



with 



$n (N, Z) = {N - 82) (126 -N){Z- 82) (126 - Z) , (6) 



^ 20 (N, Z) = (N- 82) 2 (126 - A^) 2 (N - 104) (7) 



for region A ||10|, and 



t\,, f (N.Z) ■ (^-) [^, l i>, i {X.Z) + ^ i i>-, i (X Z)+ ll2 .Y 12 (.V.z)] (8) 



with 



$21 (N, Z) = (N- 126) 2 (184 - X) 2 {Z - 82)(126 - Z) , (9) 

$ai (AT, Z) = (N- 126) 3 (184 - JV) 3 (Z - 82)(126 - Z) , (10) 

Xi 2 (JV, Z) = (N- 126)(184 -X)(X - 155)(Z - 82) 2 (126 - Z) 2 (Z - 104) (11) 

for region B | |10| |. 

In eqs. (|)-(D and (§) A = X + Z and T = \T Z \ = ±|JV - Z| 0. The respective values 
of (A^o,Zo,A ) in regions A and B are (82, 82, 164) and (126, 82, 208). The coefficients 
multiplying the functions of N and Z are adjustable parameters determined by least squares 
adjustment to the data, separately for region B [10] and for region A [10,12]. Their values 
are given in the table. 

We discuss extrapolatability first for region B. The experimental data used in the adjust- 



ment [IIJ included 211 masses. Fig. 2 shows the deviations from the data of the predictions 
of eq. (|l|) for the 56 presently known new masses measured after the adjustment flQfl . The 
respective average and rms deviations are 53 and 236 keV, as compared to 2 and 126 keV 
in ref. [pL0|1 . The corresponding deviations of Q a values are —3 and 220 keV, as compared to 



-6 and 162 keV. 

The deviations of the seven N = 126-128 nuclei, denoted by empty circles in the figure, 
increase when Z increases along the common boundary of regions A and B away from the 
data. They are related to the increasing discontinuity of the extrapolated mass surface 
along the common boundary N = 126 of regions A and B away from the data, when the 
two regions are adjusted separately PJTH|]. The deviations of the remaining 49 nuclei with N 
> 129, which do not follow the N = 126 boundary but extend into the interior of the shell 
region, seem more random and they are smaller, with respective average and rms deviations 
of —1 and 155 keV. 

The above deviations are as a rule about twice smaller than those of several recent mass 
models. This is presumably due mainly to the inclusion in eq. ([[]) of the particle- hole {p — h) 



•tin the as yet unknown odd-odd N = Z translead nuclei the ground state (g.s.) is expected to 
have T = \T Z \ + 1 and seniority zero, whereas eq. (||) with T = \T Z \ gives the energy of a low excited 
seniority two state [O] 



symmetric configuration interaction terms Ed e f, eq. (§). A more detailed discussion is given 



in ref. 11 



TABLE I. Values of the coefficients of eq. (1) determined by adjustment to the data. 



Region B 



Coefficient 


Value (keV) 


a 


-2.3859605 x 10 6 


Pi 


-1.496441 x 10 4 


fa 


-3.3866255 x 10 4 


71 


3.022233 x 10 1 


12 


2.811903 x 10 1 


73 


-3.6159266 x 10 2 


Oi 


8.16 x 10 2 


e 2 


1.007 x 10 3 


/! 


-2.121 x 10 2 


a c 


8.111517 x 10 5 


P c 


2.0282913 x 10 4 


1 C 


1.0930065 x 10 2 


¥>21 


-9.87874 x 10~ 5 


9^31 


3.13824 x 10~ 8 


Xl2 


-1.428529 x 10~ 7 



Region A (lCJ.[l2] 
Coefficient Value (keV) 



a 


-1.987628 x 10 6 


P 


-2.4773664 x 10 4 


7 


-8.51085 x 10 1 


£ 


4.585496 x 10 2 


e 


1.2183 x 10 3 


K 


2.1937 x 10 3 


a c 


7.968418 x 10 5 


P C 


2.032906 x 10 4 


1 C 


9.819137 x 10 1 


¥>n 


-4.794 x 10~ 2 


^20 


9.095 x 10- 4 



The situation in region A is less simple. The experimental data used in the adjustment 
TIJ included 29 masses and 62 Q a values connecting unknown masses. The respective 
average and rms predicted deviations for the presently known 121 new masses which became 
available after the adjustments increase drastically to —807 and 1008 keV, as compared to 
—29 and 146 keV in ref. fCDf . For the 31 new Q a values, though, the respective deviations 
are only 40 and 89 keV, as compared to 5 and 103 keV. 
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FIG. 2. Deviations of the mass predictions eq. (1) from the data for the 56 new masses in region 
B measured after the original adjustments were made. 



In order to restore to the new mass predictions the same quality as the old predictions 
had, while at the same time retaining the high quality of Q a predictions, we made |12| a 
least-squares adjustment of eq. ([I]) to all the 150 known masses, with only four adjustable 
parameters a,e,Q and k (eq. (|3|)) which largely cancel in Q a , while the other seven coeffi- 
cients were held fixed on their old values |TIJ. These are the values given in the table. 

The readjusted values of a, and k are higher and that of e is lower than in ref. 



10 



indicating smaller overall binding, smaller symmetry energy coefficient, and increased pairing 
energies in proton-rich nuclei away from stability. 

Fig. 3 shows the deviations from the data of the predictions of the readjusted eq. (|T|) for 
all the 150 experimentally known masses. The respective average and rms deviations are 2 
and 246 keV. The corresponding deviations of the predicted Q a values are 2 and 99 keV. 

A more detailed discussion is given in ref. ||12||. 
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FIG. 3. Deviations of the mass predictions eq. (1) from the data for all the 150 presently known 
masses in region A. 



III. APPLICATIONS TO SHE 



Panel a of fig. 4 shows the chain of measured a-decay energies starting from 293 118 [|TJ, 
and the values predicted for them by eq. ([!]) when the data are interpreted as g.s. transitions 
of the assigned nuclei. The figure shows as well the predictions |14| which motivated the 
search undertaken in Jl|. The respective average and rms deviations of the predicted values 
from the data are —197 and 308 keV for eq. ([T]) and —154 and 357 keV for ref. Hj. The rms 
deviation of eq. (|l|) is consistent with the deviations of the new nuclei in fig. 2 considered 
above, but the average deviation is too negative. The respective average and rms deviations 
of the predictions of eq. ([I]) from those of ref. [[14] 



are —43 and 369 keV. 



The variation of the predicted values of eq. ([!]) along the chain is smoother than that 
of the data, where there are kinks at Z = 112 and 116 presumably representing submagic 
numbers or other structure effects. In the SSME such effects are assumed [|,0] to have been 



smoothed out by configuration interaction, represented by the terms E^f. The SSME is 
inadequate to describe non-smoothed abrupt local changes associated with subshell structure 
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FIG. 4. Experimental [1] and predicted Q a values of the 293 118 decay chain, (a) Predictions of 
eq. (1) and of Smolahczuk [14]. (b) Predictions of Cwiok et al. [15], Moller et al. [16] and Aboussir 
et al. [17]. 



On the other hand, the microscopic energies calculated in ref. |14[ are basically sums of 
(bunched minus unbunched) single nucleon energies, and as such have (magic and) submagic 
gap effects built in. The corresponding predicted line in fig. 4 has kinks at Z = 110 and 
116, corresponding to predicted submagic numbers Z = 108 and 116. 

Most of the smoothing effect of configuration interaction is missing in macroscopic- 
microscopic Strutinsky type and in self-consistent mean field calculations. The included 
T = 1, J = pairing correlations seem not to be enough. This might result in calculated 
submagic gaps and associated kinks which are too large compared to the data. Panel b of 
fig. 4 shows such predicted large kinks |T7|JT^JT5[ . 

As a second application we address the a-decay chain observed in ref. [[Tj|, which is 
considered a good candidate for originating from 289 114 and its sequential decay to 285 112 
and 281 110. The respective average and rms deviations of the predictions of eq. ([!]) from the 
measured energies are 847 and 905 keV, which considerably exceed the deviations expected 
from fig. 2 for g.s. transitions. If the above assignments are confirmed, the large deviations 
might indicate that the decay chain does not go through levels in the vicinity of the g.s. 

It might also be worthwhile mentioning, that for the conceivable parents 288 112 or 291 113 
which can be obtained from the compound nucleus 292 114 by respective la or lp evaporation, 



the corresponding average and rms deviations of eq. (JH) from the measured energies are —181 
and 366 keV and —242 and 417 keV, which are more than twice smaller than for the parent 
289 114. 

In conclusion we address the bearing of the Z = 118 results on the question of magicity 
of Z = 126 versus 114. The energy eq. flip comprises a part E pair + Ec ou i, eqs. (fj) and 
(H), which is a quadratic parity-dependent function of the valence particle numbers, and a 
part Edefi eq. (|8|), which is a p — h symmetric function of both particle and hole numbers. 
For different assumed upper shell boundaries the hole numbers are different and E^f is a 
different function of N and Z. The superior agreement of eq. flip with the data in fig. 4, as 
compared to the deterioration of the SSME |J near Z = 114 mentioned in the introduction, 
demonstrates the superiority of proton-hole numbers defined with respect to Z = 126 as a 
proton magic number rather than Z = 114. This conclusion has also been arrived at in ref. 
Qj on the basis of B(E2) systematics. It is important to emphasize, though, that the above 
rather suggestive results are not a proof of superior magicity of Z = 126 as compared to the 
recently predicted Z = 120 or 124, because no comparative mass studies of this kind were 
made. 

We thank Stelian Gelberg and Dietmar Kolb for help with the calculations. 
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